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ABSTRACT 

Many toxin-antitoxin (TA) loci are known to strongly 
repress their own transcription. This auto-inhibition 
is often called 'conditional cooperativity' as it relies 
on cooperative binding of TA complexes to operator 
DNA that occurs only when toxins are in a proper 
stoichiometric relationship with antitoxins. There 
has recently been an explosion of interest in TA 
systems due to their role in bacterial persistence, 
however the role of conditional cooperativity is still 
unclear. We reveal the biological function of condi- 
tional cooperativity by constructing a mathematical 
model of the well studied TA system, relBE of 
Escherichia coli. We show that the model with the 
in vivo and in vitro established parameters repro- 
duces experimentally observed response to nutri- 
tional stress. We further demonstrate that 
conditional cooperativity stabilizes the level of anti- 
toxin in rapidly growing cells such that random in- 
duction of relBE is minimized. At the same time 
it enables quick removal of free toxin when the 
starvation is terminated. 

INTRODUCTION 

Toxin-antitoxin (TA) loci are present in many bacteria 
and archaea (1). Toxin normally inhibits cell growth, 
whereas antitoxin neutralizes the activity of toxin by 
forming a tight TA complex. 

One of the known functions of TA loci is to respond to 
nutritional stress, namely, toxins are activated upon nutri- 
tional starvation and slow down the rate of translation (2). 
Another significant feature of TA loci is that they contrib- 
ute to persister cell formation in growing bacterial cultures 



(3-5). Persisters are cells that have entered a slow-growing 
or dormant state in which the cells are tolerant to envir- 
onmental insults such as antibiotics; thus persisters 
are multidrug tolerant and therefore pose a medical 
problem. Especially, the recent experiments by 
Maisonneuve et al. (5) demonstrated that successive 
deletion of 10 mRNase-encoding TA loci of Escherichia 
coli progressively reduced the level of persisters. TA loci 
have multiple complex levels of regulation involving both 
positive and negative feedbacks and sequestration through 
binding. The importance and role of these regulations is 
still an open question. 

The relBE locus of E. coli is one of the best studied TA 
model systems. The relBE locus encodes for antitoxin 
RelB and toxin RelE. RelE is an mRNase that cleaves 
mRNA positioned in the ribosomal A site (6), including 
its own mRNA, while RelB inactivates RelE by forming a 
tight complex with it (7). RelB is a metabolically unstable 
protein whereas RelE is stable (2). However, RelB 
is translated at a higher rate than RelE, and in exponen- 
tially growing cells the abundant RelB molecules 
{[RelB] ss 10[RelE] (8)} will quench RelE activity 
completely. 

It has been shown that RelB and the RelB-RelE complex 
autoregulate relBE transcription in a complex way (9) 
(Figure 1): if only RelB is present then a RelB dimer 
(RelB 2 ) will repress relBE transcription. When RelE is 
present at a concentration such that [RelB 2 ] > [RelE] then 
a RelB 2 RelE complex binds strongly and cooperatively 
to the relBE promoter and represses transcription (9). 
In contrast if RelE increases such that [RelE] > [RelB 2 ], 
then the excess RelE molecules will destabilize the 
RelB 2 RelE-operator complex and thereby induce strong 
transcription from the relBE promoter (8). This sensitivity 
to the proper ratio between RelBE proteins is called condi- 
tional cooperativity (9, 10). 
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Figure 1. (A) Model description: antitoxin RelB (B) and toxin RelE (E) are encoded on the same mRNA (m), with only 1/100 of the ribosome that 
translate relB continues to translate relE. RelE cleaves mRNAs when it is free. RelB forms dimer RelB 2 when it is free, and RelB 2 and RelE can form 
two kinds of complexes: RelB 2 RelE and RelB 2 RelE 2 . For simplicity, we use RelB 2 as a unit for RelB and do not consider RelB monomers explicitly. 
Translation rates for RelB 2 is trans B = 15/min, and for RelE is trans E = 0.3/min. In our simulation there will be 44 nM RelE and 200 11M of RelB 2 in 
total in the non-starved condition. A list of parameters and references used are given in the 'Materials and Methods' section. (B) Visualization of 
conditional cooperativity due to the formation of RelB 2 RelE 2 , that does not repress the promoter. 



Conditional cooperativity has been found in all 
plasmids and chromosome-encoded TA loci investigated 
so far, including relBE of E. coli (8), vapBC of Salmonella 
enterica (11), phdjdoc of plasmid PI (12) and ccdA/ccdB of 
plasmid F (13). These systems belong to evolutionary in- 
dependent families (14) and function by different molecu- 
lar mechanisms. Thus conditional cooperativity must have 
an essential role in the biology of the TA genes. This 
role is difficult to analyse because of the complex 
interweaving of the components that control TA operon 
transcription. 



The relBE locus is also one of the first systems whose 
response to nutritional stress has been studied in detail (2). 
It has been found that exponentially growing wild type 
E. coli cells that are starved for amino acids reduce their 
global translation rate to a new steady state level of ~5% 
of that of the non-starved level. This dramatic change 
occurs within 20 min after starvation. Interestingly, 
deletion of relBE instead results in a post-starvation 
level of translation of ~10%. This 2-fold increase in trans- 
lation is consistent with RelE being a global inhibitor of 
translation. 
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In an attempt to understand the biological role of con- 
ditional cooperativity and to analyse the TA operon tran- 
scription in general we present a mathematical model of 
relBE operon activity that takes into account the known 
features. We subsequently investigate the role of condi- 
tional cooperativity, and show that it provides a mechan- 
ism to stabilize the level of antitoxin in rapidly growing 
cells such that random induction of relBE is minimized. 
Another important prediction is that conditional 
cooperativity enables quick recovery from the RelE- 
mediated reduction of translation when the starvation is 
terminated. 



MATERIALS AND METHODS 

In our mathematical model of the relBE system, the tran- 
scription of the relBE operon (production of mRNA) and 
translation of relBE mRNA (proteins production) are 
taken into account as two separate processes. RelB form 
tight dimers, hence we can assume that RelB is always 
present as a dimer. In addition, RelB 2 and RelE can 
form complexes in the two stoichiometric forms 
RelB 2 RelE and RelB 2 RelE 2 (32). 

A key goal of our model is to investigate the complex 
autoregulation feedback of relBE promoter activity. RelB 2 
represses transcription. With moderate amount of RelE 
added, this repression become stronger because the 
RelB2RelE has larger affinity to the operator. When 
RelE is further increased, however, the promoter is 
de-repressed because RelB 2 RelE 2 does not bind to the 
operator. 

RelE works as an mRNase only in free form, by 
degrading the relBE mRNA as well as all other mRNAs 
in the cell. In our modeling, only the cleavage action of 
RelE on relBE mRNA is directly taken into account. 

Free RelB 2 has a very short half-life since it is actively 
degraded by Lon protease. RelB half-life was measured to 
be ~3min (9), whereas RelE is stable and its half-life is 
equal to the average E. coli doubling time (~30 min in 
normal growing condition). We do not model the cell 
division explicitly, it only enters implicitly into the RelE 
half-life. 

Based on the fact that roughly 10 times more RelB 
monomers than RelE monomers are present during 
steady state cell growth in spite of the 10-fold difference 
in the half-life, we assume that the translation rate of RelB 
(monomers) is 100-fold higher than that of RelE. 

RelB 2 in complex with RelE is known to be relatively 
stable (8). However, for RelE to become active, RelB in 
complex with RelE must be degraded at some rate. We 
assume RelB 2 in RelB 2 RelE and RelB 2 RelB 2 complexes 
has a half-life roughly 3-4 to fold longer than the free 
RelB 2 (~12 min). 

The half-life of relBE mRNA is not known, and in our 
model it is set to 5 min. This is on the long side of typical 
E. coli mRNA half-life (29), which helps to keep the 
maximal promoter activity and the translation rate 
within biologically plausible values while having enough 
proteins. 



The behavior of the system is investigated throughout 
three different phases. 

A first phase, from time 0 to time 200 min in the plots, is 
what we call the 'non-starved state', where all the param- 
eters used in the simulation refer to the exponential 
growth phase of the cell. 

At time 200 min we switch to the nutritional stress phase 
(amino acid starvation). Within the framework of our 
model this means a sudden decrease in translation rate 
for both RelB and RelE by 10-fold, based on the obser- 
vation that a relBE~ deletion strain shows a reduction of 
translation to a post-starvation level of 10%. Because the 
dilution by cell division does not happen in this phase, the 
half-life of t e is changed to 24 h, which is much longer 
than the examined amino acid starvation duration (5h). 
We investigate whether enough free RelE can be released 
upon nutritional stress, since this circumstance could 
explain the 2-fold difference of the translation rate in 
starvation between wild type and relBE~ strain (2). In 
addition, starvation is known to significantly increase 
Lon activity (17), thus during starvation the half-life of 
the RelB 2 is reduced by a factor of 8, both in the free 
form and in complex. 

At time 500 min, we then switch-back to the non- 
starved set of parameters corresponding to refeeding of 
the cells with amino acids and monitor the recovery of 
the system during the switch-back to exponential growth. 

Note that the change of the parameters at the switching 
of the states has been done instantaneously. We discuss 
the effect of a time delay in the parameter change later. 

We used the Gillespie algorithm (15) and performed 
stochastic simulations (simulation procedures given 
below). relBE mRNA, RelB 2 , RelE, RelB 2 RelE and 
RelB 2 RelE 2 are the molecular players in the simulations. 
The concentrations are converted to number of molecules 
so that 1 nM corresponds to 1 molecule in a cell, which is a 
typical estimate based on the size of E. coli. Each chemical 
reaction event happens at a random time and it is chosen 
according to the reaction rates in Table 1. The possible 
reaction events are listed in Table 2. 

In the results in the main text, we consider the presence 
of four chromosomes in each cell, that means four relBE 
promoters, which is an average number of chromosomes 
for exponentially growing E. coli cells. We also tried 
the one chromosome case, which exhibit increased noise, 
but the average trajectories remains similar. Each chromo- 
some has two operators: each of them can be bound either 
by a RelB 2 RelE complex or RelB. Since cooperativity in 
the interaction between RelB 2 and the relBE operator is 
not proven (9) we assume that only one RelB 2 can bind to 
the promoter at a given time, while either one or two 
RelB 2 RelE can be bound to the operator (9). When 
the promoter is free, it shows maximal promoter 
activity, and ao/4 relBE mRNA {per promoter} per 
minute are produced. When either one RelB 2 is bound 
to the promoter or two RelB 2 RelE bind cooperatively, 
the promoter is repressed and no relBE mRNAs are pro- 
duced. In the present simulation, we did not consider the 
repression by one RelB 2 RelE because experimentally 
the Hill coefficient close to 2 in repression is observed 
in the wide range of RelB 2 RelE concentration (9). 
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Table 1. Set of parameters 



Symbol and meaning 


Description 
Value 


Units 


Reference 




a 0 total promoter activity 

[#2]ss steady state total concentration of RelB dimers 


154.665 
200 


nM min~' 
nM 


See text 
See text c.f. 


(8) 


[£] ss steady state total concentration of RelE 


44 


nM 


See text c.f. 


(8) 


t ™RN A tinlf lifp 


7 1 




kJCC LCXL C.l. 




t B RelB half-life 


4.3 


min 


See text c.f. 


(8) 


x E RelE half-life 


43 


min 


See text 




x c RelB 2 half-life in complexes 


17 


min 


See text c.f. 


(8) 


n H Hill's coefficient 


2.3 




(8) 




trans B RelB translation rate 


15 


min -1 


See text 




trans E RelE translation rate trans B /50 


0.3 


min~' 


See text 




/rbind binding on-rate 4n DajV^n 


3.S 


. _ i 
mm 


(28) 




K d B2E dissociation constant for B2E complexes formation 


0.3 


nM 


(9) 




^§=K lt B2E 

KjB2E2 dissociation constant for B2E2 complexes formation 


0.3 


nM 






k u B2E dissociation rate for B2E 


1.14 


nM 






k B2F = fei.- i x KjB2F 










k u B2E2 dissociation rate for B2E2 


1.14 


nM 






k u B2E2 = k bmll x K d B2E2 










K m dissociation constant for B binding to DNA 


10 


nM 


See text c.f. 


(16) 


K D2 dissociation constant for second B2E bound to DNA 


0.04 


nM 


See text c.f. 


(16) 


K m dissociation constant for first B2E binding to DNA 


30 


nM 


See text c.f. 


(16) 


k c cleavage rate 


2.0 


nM -1 mirT 1 


See text c.f. 


(6) 



Table 2. Events and rates in the simulation 



Event 



Rate 



mRNA transcription 

RelB dimers translation 

RelE translation 

relBE mRNA degradation 

RelB dimer degradation by Lon 

RelE degradation due to cell dilution 

RelB 2 RelE formation 

RelB 2 RelE 2 formation 

RelB 2 RelE dissociation 

RelB 2 RelE 2 dissociation 

degradation of RelB 2 RelE due to cell dilution 
degradation of RelB 2 RelE 2 due to cell dilution 
degradation of RelB 2 in RelB 2 RelE complex 
degradation of RelB 2 in RelB 2 RelE 2 complex 
binding of RelB 2 to operator 
binding of RelB 2 RelE to operator 
unbinding of RelB 2 from operator 
unbinding of one out of two RelB 2 RelE from operator 
unbinding of single RelB 2 RelE from operator 
cleavage of relBE mRNA 

stripping of RelB 2 RelE complex bound to operator from it by RelE 



^f* no. of operators with two free binding sites 
[mRNA] * trans B 
[mRNA] * trans E 

[mRNA] 

T,„ 

[RelBd 

[«m 

fcbtad * [RelB 2 ] * [RelE] 
/fbind * [ReIB 2 Re\E] * [RelE] 
k u B2E * [RelB 2 RelE] 
k u B2E2 * [RelB 2 RelE 2 ] 

{RclB 2 RclEj 

[RelB 2 RelE 2 ] 

T-E 

[ReBzBetEj 
[RelB 2 RelE 2 ] 

^bind * [RelB 2 ] * no. of operators with two free binding sites 

^bind * [RelB 2 ] * no. of operators with at least one free binding site 

(K DI * A'bind) * no. of operators with 1 RelB 2 bound 

(K D2 * /fbind) * [RelB 2 ] * no. of operators with two RelB 2 RelE bound 

(K m * Ar b i nd ) * [RelB 2 ] * no. of operators with a single RelB 2 RelE bound 

k c * [RelE] 

^bind * [RelE] * no. of operators with at least on RelB 2 RelE 
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Algorithm of the Gillespie simulation 

We use the Gillespie algorithm in our simulation, where 
the chemical reactions are treated as discrete stochastic 
events that happen at given rates with time interval 
between events drawn from the exponential distribution 
(15). 

The state of the system at each time step is defined by 
the concentrations of the five molecular players, i.e. RelBE 
mRNA (m), free RelE (E f ), free RelB 2 (B t ), RelB 2 RelE 
(B 2 E) and RelB 2 RelE 2 (B 2 E 2 ). The total copy number of 
RelB 2 (B T ) and RelE (E T ) are given by 



The last four reactions are what we call 'stripping' (18,19), 
where RelE forming complex with RelB 2 RelE bound on 
the operator and removing it, and the reverse reaction of 
the stripping. 

Each run is from time 0 to time 600. At time / = 200 
min the values of x B , t c , t e , trans B and trans E are changed 
from values that mimic fast growth conditions to values 
typical of amino acid starvation. At time / = 500 the same 
parameters are changed back to the fast growth value. 

All the presented results refer to the averages of the 
concentrations over a sample of 1000 simulations, unless 
otherwise noted. 



Bf = Bf B 2 E -)- B 2 E 2 , 
Ef = Ef -\- B 2 E -\- 1B 2 E 2 . 

The chemical reaction with a rate specified by the param- 
eters showed in Table 1 results in a change in the number 
of molecules as follows. 



ao/4 • no. of unrepressed promoters 

m ► m + 1 

m > m — 1 

trans B m+B 2 EK DB 2E -h 

B { ► Bf + 1 

Bf^Bf+l 
Ef Bf + 2 

Ef ► Ef + 1 

Ef^Ef-l 

B 2 E * r£A+i?2£rJW ^ B 2 E+ 1 
B 2 E B2EK ™ ME ^ B 2 E-\ 
B 2 E 2 — - * B 2 E 2 -)- 1 

D j-, B2E2-K D B2E2-kh+B2E 2 /i: c 
£> 2 £L 2 ' £> 2 h 2 — 1 

In addition, RelB 2 and RelB 2 RelE can bind to the four 
operator sites. The binding happens to each operators in- 
dependently, and the binding rates are given as follows 
(O expresses the operator): 



O + Bf'^^O-Bf 

O.B t KMOB!) > 0 + Bf 

O + B 2 E kh 0 BlE ~, o ■ B 2 E 

O . B 2 E W ' A ' (0,Ji2£) > O + B 2 E 

O ■ B 2 E+B 2 E IS:i:m n i: < o ■ (B 2 E) 2 



O ■ (B 2 E) 2 K ^°< B ^\ O B 2 E+ B 2 E 

O ■ B 2 E+ Ef k "- E '<°-^.. 0 + B 2 E 2 

O + B 2 E 2 k ^ K ^ K ^°- B ^- ) o . B 2 E+Ef 

O ■ (B 2 E) 2 + Ef A>£K0 - (Ji2£ H O ■ B 2 E + B 2 E 2 

O B 2 E+ B 2 E 2 O ■ (B 2 E) 2 + Ef 



(1) 



RESULTS 

Outline of the model 

The overall regulations and feedbacks in our model are 
summarized in Figure 1A. The mechanism of the condi- 
tional cooperatively is illustrated in Figure IB. RelB and 
RelE can form two types of complexes, namely RelB 2 - 
RelE (B 2 E) and RelB 2 RelE 2 (B 2 E 2 ). B 2 and B 2 E repress 
the promoter activity of relBE operon, while B 2 E 2 does 
not. This is a scenario proposed to reproduce the observed 
conditional cooperativity (9). Both RelB 2 RelE complex (9, 
32) and RelB 2 RelE 2 complex (as heterotetramer); (32) 
have been observed in vitro. 

All the molecules are exposed to either degradation by 
proteases or dilution by cell division, where the details are 
described in the next subsection. The mRNA m can be 
also actively degraded by the free toxin E. 

Model parameters 

The parameters of the model are constrained by (i) 
Stoichiometric data showing that when [total amount of 
RelE]:[total amount of RelB monomer] is in 1:2 ratio, 
RelB 2 RelE and operator O complex, (RelB 2 RelE) 2 - O, 
is formed, while increasing the ratio of RelE further 
destabilizes the complex (9). This can be reproduced 
when the binding of RelE to RelB 2 and binding of 
RelE to RelB 2 RelE occur with similar dissociation 
constants, hence when the total amount of RelE exceeds 
the total amount of RelB 2 , RelB 2 RelE are converted 
to RelB 2 RelE 2 (See Supplementary Material A for 
detail). 

(ii) That under normal growth in rich medium the in 
vivo concentration of RelE is about one-tenth of that of 
RelB (8). The actual concentration level was estimated to 
be 550- 1100 nM for RelB while 50-100 nM for RelE (8). 
Here we choose parameters so that total RelB is ~500 nM 
and the total RelE is ~50nM in non-starved exponential 
growth condition. 

The in vivo lifetime of RelB is ~4.3 min, whereas RelE is 
metabolically stable but diluted by a rate set by cell 
division, giving it a characteristic lifetime of ~43min in 
the exponential phase. This 10-fold shorter lifetime for 
RelB than RelE, with the 10-fold higher concentration 
of RelB in the exponential growing condition mentioned 
above gives us a translation rate of RelE ~ 1 % of that of 
RelB. [Since the estimate of the RelB and RelE levels in 
vivo was difficult due to the low cellular amounts of the 
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Figure 2. Response to amino acid starvation and later recovery. The system is starved for amino acids from 200 min to 500 min. (A) Probability 
distribution P([Ef\, ') of a cell having a certain concentration [£ r ] nM of free RelE at a given moment /. (B) Time courses averaged over 1000 cells, for 
free RelE, free RelB and relBE mRNA, illustrating how the system switches between a state of high antitoxin to a state of high free toxin. (C) The 
dynamics at entrance to the starvation at the single cell level. Three examples are shown, and the total amount of free RelE is plotted as function of 
time, from time 180 to time 300. 



proteins (8), the 10-fold difference in concentration can be 
smaller in reality, which would indicate higher translation 
rate ratio (up to 10%).] (hi) That RelB 2 represses RelBE 
promoter by a factor 16, whereas (RelB 2 RelE) 2 represses 
the same promoter by a factor 800 under exponential 
growth in normal medium. These repression folds are 
about one-tenth of what was observed using lacZ fusion 
on low copy number (~10) plasmids (16). We adapted 
these repression folds mainly because they are close to 
the upper limit of the repression folds which can give 



biologically plausible promoter strength for the relBE 
promoter and translation rate with still being able to 
reproduce the observed concentrations of RelB and 
RelE in exponential growth, where promoters are fully 
repressed. 

In the next subsections, we demonstrate that the model 
with these evaluated parameter values show expected 
response to amino acid starvation and recovery when 
the starvation is over. We also investigate the robustness 
of the behavior against the parameter changes. 
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Amino acid starvation drives the switch to toxin activation 

The model is examined when switching from rich medium 
to amino acids starvation, and subsequently exposed 
to refeeding after 5h. There are three main processes 
taking place during starvation: (i) The overall translation 
rate is reduced, because of the lack of amino acids, here 
simulated by a 10-fold reduction, (ii) Cells stop dividing, 
thus dilution of RelE decreases, allowing for accumulation 
of this long lived protein, (hi) Lon activity is believed to be 
increased under starvation (17), increasing the degrad- 
ation rate of RelB. We chose to make the degradation 
8-fold stronger during amino acid starvation. The result 
presented here holds as long as the degradation of RelB in 
the RelB 2 RelE complex is large enough during the starva- 
tion, and the 8-fold is close to the minimum fold needed. 

Figure 2 shows switching from the antitoxin dominated 
state to the toxin dominated state elicited by amino acid 
starvation. Importantly, in order to make this response 
work we had to assume that RelB in complex with RelE 
must be actively degraded, at least during starvation. 
If such degradation was not included, then the transition 
to the high-RelE state cannot be achieved. Here, we 
assume that RelB is degraded 4-fold less effectively when 
in complex with RelE than when it is free {in vitro data 
supports that RelB is partly protected from degradation 
when it is in complex with RelE (8)). To demonstrate the 
necessity of this process, we compare the model with and 
without active degradation of RelB in the complex in the 
Supplementary Material B. 

From Figure 2B, we also see that there is a some time- 
lag to enter the toxin dominated state, while the time to 
exit this state is very short. This time lag is due to the com- 
bined effect of 'stripping' and conditional cooperativity. 
By 'stripping' we mean the process where a free toxin 
molecules 'invades' the RelB 2 RelE complex bound to the 
promoter inducing the complex to be released from 
operator DNA. This has been shown to occur both 
in vitro and in vivo (8). Note that it is possible to have 
conditional cooperativity without stripping, in which 
case excess RelE will form RelB 2 RelE 2 complex in the 
bulk and thus sequester RelB 2 RelE, but does not 
remove bound RelB 2 RelE actively from the operator site. 

The conditional cooperativity with stripping opens for 
an active and relatively long battle between degradation of 
RelB and a de-repression of the promoter with an 
associated rise of relBE mRNA (2,8) and hence increase 
in production of RelB. [It should be noted that in the 
experiment that relBE mRNA level was observed to rise 
~30-fold just after the amino acid starvation (2,8), 
while in the present model we observe only ~2-fold rise 
(Figure 2). The height of this peak depends strongly on the 
cleavage rate of mRNA by free RelE, k c , but just lowering 
this value delays to reach the high free toxin state at the 
starvation (details in Supplementary Material C). This dis- 
agreement can be in principle improved if we take into 
account the fact that free RelE will interact with all the 
mRNAs in the cell, thus it is quite likely that free RelE will 
be sequestered and will not cleave much of relBE mRNA 
when its concentration is very low, which should give 
bigger rise for relBE mRNA. In the present model, 



however, we do not take this effect into account because 
of large ambiguity in detailed interactions between all 
mRNAs the free RelEs.] Without stripping, this rise 
becomes much smaller than this, and hence the delay 
becomes less (details in Supplementary Material D). 
Central in this 'battle' is cleavage of mRNA by free 
RelE, since it reduces the relBE mRNA for both RelE 
and RelB, and thereby favor the long-lived toxin RelE. 
The cleavage rate k c of relBE mRNA by toxin is not 
known in vivo, and we use £ c = 2/min/nM to obtain rea- 
sonably fast rise of the toxin upon starvation. For com- 
parison, the in vitro cleavage activity per codon with 
empty ribosome A-site was estimated to be between 
0.000042 to 2.4 /min/nM depending on the codon (6). 
Considering there are 79 codons for RelB, one expect 
2-3 ribosomes at any time to translate the mRNA, and 
thus an effective cleavage rate that should not exceed 
5/min/nM. Thus our assumed value is in the high end, 
but using for example a 10 times smaller value of k c 
would delay the transition into toxin dominated state by 
hours, as shown in Supplementary Material E. 

Single cell activation of RelE is binary 

The behavior of single cells are summarized in Fig. 2A 
as the probability distribution P([E t ], t) of a cell having 
a certain concentration [E t ] nM of free toxin at a given 
moment t. We can see that the response is binary. There 
is a high peak at low toxin at the start of starvation 
(time 200), but another peak for high free toxin state 
(around [.Ef] « 45 nM) appears already at 10 min after 
starvation. The low probability to take the value in 
between (10 nM < [E t ] < 30 nM) suggest that each cell 
switches from low (<10nM) to high free toxin (>30nM) 
quickly. On cell population level (Figure 2B) in contrast it 
takes almost 40 min for free RelE to reach 30 nM. This 
contrast reflects the big variation in switching time 
between different cells. 

Figure 2C shows three representative trajectories of 
single cells entering starvation. At the beginning, free 
RelE is almost zero because RelE are bound by RelB's, 
but after some time the balance switches and the concen- 
tration of free RelE rises quickly to high level. 
The examples show variation in switching timing over 
60 min but none of them spend much time at the 
intermediate free toxin level. Overall this shows the signifi- 
cance of stochastic modeling of this type of system as 
well as the need for single cell measurement of the TA 
systems. 

Amino acid starvation is terminated at t = 500 min in 
Figure 2. Compared with the entry to the starved state, the 
recovery from the starved state is found to be extremely 
quick both on average and at the single cell level. As 
demonstrated below the fast recovery depends on condi- 
tional cooperativity. 

Conditional cooperativity primes fast exit from the toxin 
dominated state 

To clarify the role of conditional cooperativity, we 
compare the system with conditional cooperativity (the 
same system as in Figure 2: referred as 'cc') and the 
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Figure 3. Role of conditional cooperativity. (A) The time evolution of free RelE level for the system with (red) and without (blue) conditional 
cooperativity. The system is starved for amino acid from 200 to 500 min. (B) Probability distribution of free RelE in the non-starved state without 
conditional cooperativity (blue) and with conditional cooperativity (red). Free RelE takes higher value without conditional cooperativity. 
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Figure 4. Summary of the model behavior against parameter changes. For each parameters (horizontal axis), fold change of the values from our 
reference values are tried one by one. The color gradients indicate how the model deviates from the reference behavior: yellow indicates too many 
free toxins in the healthy states, green indicates too slow rise of free RelE at amino acid starvation and red indicates too slow drop of toxins after the 
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system without conditional cooperativity (referred as 'no 
cc'): We removed the conditional cooperativity by pre- 
venting RelB 2 RelE 2 complex formation completely, 
while keeping the other parts of the system unchanged. 
Therefore, regardless of the relative concentrations of 
total RelE and total RelB, the only possible complex 
they can form is RelB 2 RelE. Figure 3A illustrates that 
the conditional cooperativity is required for the fast 
recovery from the toxin dominant state under starvation 



to the non-starved state. Mechanistically this reflects that 
formation of RelB 2 RelE 2 derepresses the promoter in the 
high-i\e/i? state and thereby primes the system for 
recovery already during starvation. 

Figure 3B shows that conditional cooperativity also 
reduces the probability to have high free toxin in the 
non-starved state. In the non-starved state the excess of 
RelB buffers for an increase in free RelE (by the se- 
questration into the complex RelB 2 RelE). If RelB by 
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fluctuations becomes low, the conditional cooperativity 
provides a negative feedback that secures additional 
RelB. This reduces free RelE by complex formation 
and thus the concentration of free toxin is kept low. 

Robustness of the observed behaviors against the 
parameter changes 

Finally, we study the robustness of the observed behaviors 
against parameter changes. This is not only to see how 
general our conclusion is regarding the choice of the par- 
ameter values, but also to indirectly test how the feedbacks 
that are not modeled in the present framework would 
affect the behavior. Especially, growth-rate dependences 
(24-27) are observed in cell physiology through for ex- 
ample the partitioning of the ribosomes (which affects 
the translation rate) or the RNAP availability (which 
affects the transcription rate). Such dependences will 
affect parameter values upon amino acid starvation. The 
robustness test gives idea about how significant such a 
feedback can be in the studied behavior. 

In Figure 4, we summarize the robustness of the 
observed behaviors against parameter changes. We 
change parameters (or ratio of parameters) one by one 
by 2"-fold, with — 3<«<3, and we check whether the 
model is working with that parameter value based on 
the following criteria: (i) The free RelE level in the 
healthy state is <lnM. (ii) The free RelE level reaches 
>10nM within 20 min after the start of the amino acid 
starvation, (hi) The free RelE level drops < 1 nM within 
5 min after starvation is stopped. 

The robustness analyses shows that the condition (ii) is 
hardest to satisfy. Most parameters are split in two sub- 
regions, i.e. the model is very sensitive to the change in one 
direction (increase or decrease) but insensitive to change in 
another direction. This is because our reference parameter 
is at borderline, i.e. just fast enough, to satisfy the criterion 
(ii). Even at this borderline we had to introduce the faster 
degradation of RelB's in complexes during starvation, to 
achieve fast increase in free RelE. This again support the 
necessity of the fast degradation of the antitoxins at amino 
acid starvation. 

The only case where condition (i) tends to be violated is 
when the dissociation constant K d B2E for RelB 2 RelE for- 
mation is very large, hence RelE are not tightly seques- 
tered in the complex. Increasing both K d B 2 E and K d B2E2 
has similar effect, but it affects stronger on the recovery 
from the free high-toxin state (iii), by freeing up toxin 
easier at the transient state. 

In addition, the third criterion (iii) is violated when the 
unstable antitoxin RelB is not produced high enough, 
which happens when the ratio between the translation 
rate of RelB and RelE (trans B /trans E ) is too small or the 
repression by RelB 2 RelE (characterized by K D2 ) is too 
tight. When K C ,B 2 E/K C ,B2E2 is too small (1/4-fold or 
less), RelB 2 RelE 2 is not formed as much and the condi- 
tional cooperativity becomes ineffective, which also makes 
the recovery slow. 

The robustness of the transcription rate oc 0 and the 
translation rate (trans) is of particular interest in the 
context of the growth-rate dependent feedback. We see 



that the model behavior is robust 1/8- to 2-fold change 
of these parameters. Namely, even if there is a feedback 
from the growth rate to these parameters, the model 
behavior will not be altered as long as the change is 
within this fold. Especially, the transcription rate ot 0 is 
expected to decrease upon slower growth (24), which 
is the direction where the model behavior is robust. 

We also investigated how the time scale of the param- 
eter changes (when the conditions are shifted to or from 
starvation) affect the kinetics of the transitions between 
two states. The relevant parameters changed are the RelB 
lifetime and the translation rates. Overall, we find that the 
time scale of parameter changes is rate-limiting for the 
entry into the starved state, i.e. the kinetics of the transi- 
tions is fully determined by how rapidly we change these 
parameters (Supplementary Material E). This is again 
because our reference parameters are at the borderline 
to reproduce the fast entry. 

The recovery phase is, however, less sensitive to the time 
scale of the parameter changes. The translation rate has 
the biggest effect, but as soon as the translation rate in- 
creases by some amount, the RelB accumulate enough and 
free RelEs are repressed (Supplementary Material E). 
Furthermore, the main result regarding the conditional 
cooperative is robust: the recovery is always much faster 
for the model with conditional cooperativity than without 
conditional cooperativity (Supplementary Material F). 



DISCUSSION 

We constructed a mathematical model of how relBE is 
regulated, with the main focus centered around the condi- 
tional cooperativity in the autoregulation of relBE operon. 
With our model we inteded to capture the available 
experimental data, test the known and estimate the 
values of unknown parameters and investigate the 
systems dynamics when the cells are shifted between 
non-starved and starved states. 

Our current modeling framework highlights several 
interesting features: 

(1) A fast entry to the high-toxin state can only be 
realized if antitoxin is degraded both when it is free 
and in complex with toxin during amino acid 
starvation. 

(2) The transition from the antitoxin dominated state to 
the toxin dominated state upon sudden amino acid 
starvation is not graded but binary at the level of 
single cells. 

(3) When amino acid starvation is terminated condi- 
tional cooperativity mediates fast recovery from the 
toxin dominated to the growing state. 

(4) Conditional cooperativity also reduces the occasional 
occurrence of high free toxin state in the non-starved 
condition. 

In the following we discuss these four features. Active 
degradation of antitoxin in the TA complex during the 
starvation [Feature (1) above] should be closely coupled 
to Lon protease activation during amino acid starvation 
(17). More detailed understanding of Lon activity during 
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the starvation as well as fluctuations of Lon activity in 
growing cell will provide more insights about toxin acti- 
vation through this pathway, which can be included in the 
future development of the model. 

The switch-like activation of the toxin at the single cell 
level [Feature (2) above] comes from the positive feedback; 
when free toxin starts to increase, the cleavage-rate of TA 
mRNA increases, that in turn reduces the amount of anti- 
toxin and thereby forces more free toxin to accumulate 
(see also Figure 1). Positive feedback facilitated switches 
have been seen in many other systems that require decision 
making, for example, the lysis/lysogeny switch in temper- 
ate phage (20), sugar utilization in bacteria (21) and cell 
differentiation (22). In the present case, the decision that a 
cell needs to make is whether it should translate or shut 
down translation; we believe that it makes biological sense 
that the cells do not waste time between these two states. 

Conditional cooperativity [Feature (3)] facilitates the 
switch to work in a robust manner, favoring the antitoxin 
dominated state by making the switching back dynamics 
fast and by reducing the probability to randomly switch to 
toxin dominated state without stress [Feature (4)]. 

Even though the present modeling relies on parameters 
measured for re/BE systems, we believe that our model is 
relevant for the mRNase TA systems in general, where the 
basic framework of the regulations is believed to be similar 
to the one in RelBE. Especially, the observed switching — 
between low- and high-toxin states upon starvation and re- 
covery — is interesting in relation to persisters. A current 
view on the persister mechanism suggest that by toxin 
fluctuations some cells happen to end in a toxin dominated 
state [cf. (4,23,24)], while the present study shows that if 
the TA system has the conditional cooperativity such fluc- 
tuations will be strongly suppressed (Figure 3B). The 
extent to which the system is subject to fluctuations, i.e., 
how often a cell can be in high-toxin state by chance 
without amino acid starvation, is an important quantity 
to study in the future in relation to the persister formation. 
It is in general an interesting theoretical and experimental 
problem to understand the role of TA system in persister 
cells formation in the light conditional cooperativity. 

In order to generalize the present model, especially to 
the persister system, it should be noted that the several 
known feedbacks are not taken into account in the 
present model. As mentioned earlier, a number of param- 
eters are growth-rate dependent (24-27) due to e.g. RNAP 
availability, stringent response via ppGpp, DksA regula- 
tion, etc. In the present level of modeling, these factors 
would mainly affect the transcription rate. Within the 
interest of the present work, the conclusions remain quali- 
tatively the same for changing the transcription rate by 
1/8- to 2-fold (Figure 4). However, the feedback where 
high-toxin state imposes slow growth and slow transcrip- 
tion which in turn favors high-toxin state opens for an 
interesting direction as it can have a strong effect on the 
stability of the persister state (24). 

Finally, we propose several possible experiments based 
on the present results. 

The direct test of the predictions of the model about the 
conditional cooperatively would be to construct a mutant 
that does not form RelB 2 RelE 2 and yet keep other 



properties of RelB and RelE, and then compare the 
dynamics with the wild type scenario. This, however, 
requires the detailed structural knowledge of the proteins. 

The observation of the dynamics of the entrance to and 
recovery from the starved state at the single cell level will 
give a lot of information about the system. Especially, the 
binary response [Step (2)] should be confirmed experimen- 
tally, by for example visualizing the RelE level or relBE 
mRNA level in each cell. 

The low copy number of RelE (~50nM) or relBE 
mRNA (~a few molecules) makes it challenging to 
monitor these molecules using usual fluorescent micros- 
copy, to overcome this challenge superresolution micros- 
copy (e.g STORM) has been successfully used to monitor 
the single cell/single molecule dynamics in bacterial cells 
(31). The dynamics of the recovery is also interesting to 
observe, especially the expected short recovery time due to 
the conditional cooperativity. 

Another interesting experiment is to study the depend- 
ence of recovery dynamics on the duration of the amino 
acid starvation. The duration of typical experiments is 
~5 h, and the toxin RelE is expected to be stable on that 
time scale. This is consistent with (7), where a pulse of 
RelE found to sufficient to prevent cell division for 
similar time scale. The limit of the stability of RelE in 
non-dividing cells, though, is not known. This factor, ac- 
cording to our model, can have a strong effect on the 
recovery behavior from high-toxin state after very long 
lasting amino acid starvation. 

SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online: 
Supplementary Data A-F and Supplementary 
References [33,34]. 
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